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ABSTRACT 


Wind fields were measured with the ground-based NASA/MSFC Ildar 
and with the NASA RB-57 instrumented aircraft. The remotely sensed 
winds are compared with the in situ aircraft measurements. The mean 
wind fields, the turbulence intensities, and the turbulence spectra 
determined from measurements by both systems are In very good agree- 
ment. Turbulence instensities and spectra were calculated from the 
fluctuations with time in the radial wind speed component. It should 
be noted that time histories of the radial wind represent values 
spatially averaged over a 300-m volume element. The Ildar winds were 
sampled at approximately 2 times per second whereas the aircraft 
measurements were sampled at 40 times per second. 

The second moment or Doppler frequency sepctral width of the Ildar 
measurements was also compared with turbulence msities measured by 
the aircraft. These second moments could only be resolved at the very 
low altitudes (in three range bins). Turbulence intensities estimated 
from the spectral width data were an order of magnitude higher than 
that measured by the aircraft. 

An interesting boundary layer evolved during the progress of the 
experiment. The breakup of a stable boundary layer resulted in winds 
blowing in one direction above 600 m msl and in the opposite direction 
below that level. Both the aircraft and the lidar systems clearly 
identified this unusual boundary layer flow and showed the identical 
trends. The clear identification of the unusual boundary layer by 
both systems further augment the reliability of the remotely measured 
wind speeds and turbulence. 
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1.0 INTRODUCTION 


A field test was carried out to compare Ildar-measured winds and 
turbulence; with both aircraft measurements and tower array measurements. 
The instrumentation consisted of the NASA/MSFC Ildar (Bllbro and Vaughan 
1978), the NASA RB-57 Instrumented aircraft (Camp et al. 1983), and the 
NASA/MSFC Atmospheric Boundary Layer Facility eight-tower array (Frost 
and Lin 1983). The experiment called for three days of testing. 

On May 10, 1983, the Doppler Ildar was set In a conical scanning 
mode. Scans were carried out at vertical angles of 6°, 19°, 26°, and 
32°. The aircraft then flew circular flight paths at increasingly 
higher altitudes in order to approximately capture the Ildar beam as 
illustrated In Figure 1. On May 11 the Ildar was down for repairs and 
adjustments and only flights over the tower array were carried out. On 
May 12 the lidar was fixed at a 6® vertical angle and at 52° azimuthal 
from true north, see Figure 2. The aircraft then flew approach paths at 
an approximate 4° glide slope parallel to the radar beam. Eight succes- 
sive runs or approach paths were flown at approximately 5-minute Inter- 
vals. On all three days, May 10-12, the tower array was operated. 

The aircraft made several passes directly over the towers to provide 
data for comparing flight measurements with ground-based tower measure- 
ments. These data, however, are not analyzed to any extent In this 
report. 

The emphasis of the study was to compare Doppler-1 Idar-measured 
winds and turbulence with aircraft measurements. Primarily the study 
was to compare aircraft-measured turbulence intensities with the lidar 
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second moment or spectre! width data, Unfortunately, this aspect of the 
study was not particularly successful In view of the fact that only 
three range bins (Range Bin 9, 10, and 11) had high enough signal -to- 
nolse ratio for the second moments to be successfully computed. 

Secondly, the values computed in the range from 1.26 to 2.51 m/s, 
which is a factor of ten larger than those values measured either with 
the aircraft or with the tower array. 

The field study was successful, however, in that It: (1) provided 
a unique set of data for comparing mean wind speed values; (2) revealed 
that turbulence intensities computed from the Doppler-measured wind 
speed time histories (i.e., 300 m spatially averaged values) agree 
remarkably well with the point measurement from the aircraft; and 
(3) showed that turbulence spectra calculated both from the time 
histories of the lidar-measured winds and the aircraft-measured winds 
were in very good agreement. 

Finall", an extremely interesting .atmospheric boundary layer event 
evolved during the time period (16:42-17:78 Zulu) of the May 12, 1983, 
test. This event wa M cl^-'ly recorded by both the aircraft instrumen- 
tation and the lidar. Because both systems accurately recorded this 
boundary layer event, it is believed that considerable reliability in 
the Ildar mean winds is demonstrated. 

This report presents a detailed analysis of the winds measured 
during the evolution of the atmospheric boundary layer occurring on 
May 12, 1983, and emphasizes the validation of the Doppler Ildar remote 
measurements with the in situ aircraft measurements. 
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2.0 INSTRUMENTATION AND DATA 


2.1 Lidar 

A complete description of the NASA/MSFC Doppler lidar is provided 
in Bilbro and Vaughan (1978), Jeffreys and Bilbro (1975), and Lee (1982). 
The lidar is a variably pulsed C0 2 Doppler lidar. During this study, a 
2-us-pulsed lidar was used. The Doppler measures the component of the 
wind along the lidar beam, i.e., the radial wind speed component. The 
measurements are representative of the average wind speed within a 
conical trapezoid of 300 m in length and of diameter associated with 
the diverging lidar beam width. Figure 2 illustrates the lidar beam and 
shows the location of each individual range bin for which radial wind 
speed components are measured. The figure also illustrates the position 
of the beam relative to the terrain contour cross section, 

The lidar data were received from NASA/MSFC in digitized format on 
magnetic tapes. The data format for a given tape is shown in Table 1. 
Data for roughly 189 range bins is calculated every half second, as 
illustrated in the table. Typical time histories of the data provided 
on the tape, which includes amplitude of the signal in decibels, radial 
wind velocity in meters per second (m/s), and second moment (lidar 
width) data for turbulence intensities in meters per second, are shown 
in Figure' 3 for the May 10 and 12 field tests, respectively. The lidar 
width data is recorded as an integer index. The corresponding value of 
each integer index is given in Table 2. 

Figure 4 is a plot of 150 sequential values of wind velocity from 
the May 12, 1983, data tape. The figure illustrates approximately 75 
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seconds of data. It is clear from the figure that data in Range Bins 1 
through 8 are very noisy due to ground clutter and do not provide useful 
data. Also, the figure shows that beyond approximately Range Bin 21 the 
signal -to-noise ratio becomes excessive and velocities measured above 
this altitude are not meaningful. Thus for the May 12 field test, only 
radial wind speed values from Range Bin 9 (460 m msl) to Range Bin 21 
(840 m msl) were selected for analysis. 

2 . 2 Aircraf t 

Data from the RB-57 flights consisted of 80 variables in a 60-bit 
integer format. The original raw data were sampled at 200 cycles per 
second. However, they were provided from NASA/Langley in engineering 
units at 40 samples per second. Although all the variables necessary to 
resolve the wind speed components by backing out the aircraft motion are 
available, the data from Langley provided pre-c'omputed gust velocities. 
These were used throughout the analysis. Table 3 provides a sample of 
the aircraft data. 

Data sets for eight flights along an approximate 4° glide slop* 
parallel to the lidar beam were collected. Additionally, four level 
flights perpendicular and parallel to the tower array were made and data 
provided. The tower data, however, are not analyzed to any extent in 
this report. Figure 2 shows typical flight paths relative to the lidar 
beam. Because Of unusual drift in the INS, the latitude and longitude 
measurements are questionable. Thus the exact position of the aircraft 
relative to the lidar beam in a horizontal plane is not known. Ground- 
based personnel, however, observed the aircraft to approach essentially 


4 


along the lidar beam. Moreover, the aircraft height at any instant is 
accurately measured and is, in fact, the most important value of air- 
craft position for comparing the wind speeds measured by the two systems. 

The horizontal wind measured by thtf aircraft in terms of wind 
speed, magnitude, and direction are plotted in Figure 5 for all eight 
runs. One observes that during the May 12 field test a strong inversion 
developed at approximately 600 m above the elevation of the lidar site. 
This resulted in decoupling of the wind with the wind blowing in one 
direction aloft and another direction at the surface; along some 
flight paths, the wind is observed to change direction by as much as 
180°. 

Figure 6 is a three-dimensional plot of the horizontal winds 
measured with the aircraft along each flight path and staggered in time. 
In this plot the wind vectors illustrated are values averaged over a 300 
m section along the flight path. One observes the growth of the inver- 
sion layer over the 30-minute period during which the eight flights 
were carried out. 

The temperature variation with height along the flight path was 
computed and is plotted in Figure 7. It is apparent from the data that 
although the temperature gradient aloft represents a heit flux toward 
the surface of the earth, the potential gradient is, in general, posi- 
tive, Values of the gradient and flux Richardson's number, repsectively, 
were computed based on parameters averaged along the entire flight path. 
The results were not meaningful in that the flux Richardson number was 
a Rrelatively large negative value while the gradient Richardson number 
was positive and on the order of 0.25, which represents an extremely 
stable boundary layer. 
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In order to attempt to resolve this anamoly, the velocity profile 
was broken up into multiple sections as illustrated in Figure 8. A 
linear curve was fit to each segment of the velocity profile and the 
slope (i.e., velocity gradient) for that segment computed, The values 
of the curve fit parameters are tabulated in Figure 3. The Richardson 
number was then plotted for the individual sections of the atmosphere. 
These results are plotted in Figures 9 and 10. The meaning of the 
results, however, is not clear in terms of the classical Richardson 
number for flow over flat terrain. (Note, however, that the terrain is 
not flat since Mankin Mountain (approximately 350 m msl) is directly 
under the fl ight path. ) 

Inspection of the first plot of Figure 9 is obviously not mea- 
ningful. The next two plots show unacceptably large values of Ri below 
450 m msl but neutral or slightly unstable flow aloft. The latter plots 
illustrate that the lower layer has become unstable (undoubtedly due to 
heating from the ground) and that the neutral layer is growing downward. 
The flux Richardson number given in Figure 10 shows even less decisive 
results. The reliability of the Richardson number prediction is 
obviously, at most, marginal. 

Further analysis of the data, including the results from the tower 
measurements, will be carried out under subsequent efforts to attempt to 
provide a better physical explanation of the mechanism generating the 
observed boundary layer behavior. 
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3.0 COMPARISON OF LIDAR MEASUREMENTS WITH 
AIRCRAFT MEASUREMENTS 


3.1 Mean Wind 

Comparison of the measurements of mean wind with the Ildar and with 
the aircraft system Is described in this section. The a Ire raft -measured 
wind speeds were first transformed to the time-dependent components 
alone a 6° line of sight and at 52° azimuthal true north, I.e., along 
the lidar beam. 

The aircraft-measured wind speeds were then averaged with time over 
a period corresponding to the length of time required for the aircraft 
to traverse the 300-m range bins along the flight path. Two approaches 
to carrying out this averaging technique were investigated. One was to 
assume vertical homogeneity in the flow field. The averaging process 
for the aircraft data was then carried out as illustrated in Figure 
11(a). The alternate technique was to average the wind assuming homo- 
geneity in the horizontal direction. This approach is illustrated in 
Figure 11(b). 

A third effect taken into account when comparing data from the two 
systems was to assure that the winds measured with lidar and with the 
aircraft were measured in the same time period. The run times asso- 
ciated with each flight path were therefore overlaid on the lidar- 
measured winds as illustrated in Figure 12. The lidar data are sampled 
in each bin at approximately 0.5-second intervals. The segment of the 
lidar wind speed time history associated with the time period in which 
the aircraft was passing through or parallel to that range bin was then 
averaged. 
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Figures 13 and 14 compare the lidar-measured winds averaged over 
the time period, as described above, with the aircraft-measured winds 
averaged over the corresponding 300-m section assuming vertical homo- 
geneity and horizontal homogeneity, respectively. One observes very 
good agreement between the lidar measurements and aircraft measurements 
although the data is consistently higher for the lidar measurements. 
Although the exact cause of this difference is not known, it is rea- 
sonable to assume that :ie to the unusual drift in the IN.S the aircraft 
velocity may be low because of the Schuler oscillation phenomenon. 

In general, the trends of the aircraft-measured wind most closely 
follow the lidar measurements when the assumption of vertical homo- 
geneity is made. This implies that the best agreement is achieved when 
the aircraft is at the same distance from the lidar even though it may 
be abov*’' below the lidar beam at that distance. Horizontal homo- 

geneity, of course, implies that the aircraft is making measurements at 
the same height as the lidar beam for the given range bin but may be 
further or closer to the lidar location in horizontal distance. It 
should be noted that no attempt is made to correct the velocities for 
convective effects, i.e., translation of the air pacel parallel to the 
lidar beam, nor for surface terrain contour effects. A terrain correc- 
tion may help improve the data comparison since the lidar beam passed 
directly over the top of a mountain, whereas some of the flight paths 
may have passed to one side or the other. The agreement of the data 
is believed to be sufficiently good; thus, no terrain correction was 
attempted. 
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3 , 2 Turbulenc e 

Computed turbulence intensities for the radial wind speed component 
from the aircraft measurements and the Ildar measurements are also shown 
in Figures 13 and 14. In these figures the turbulence intensity of the 
1 idar-measured wind is computed from 

* i/j, wt> - *> 2 

where Sf is the average wind speed for the period of time the aircraft 
passes through or parallel to the range bin of interest and W(t) is the 
fluctuation in wind. The summation is carried out over N time incre- 
ments' of At = 0.455 second which lapses the time interval between the 
aircraft entering and leaving the range bin. This time interval is used 
both in computing the aircraft turbulence intensity, illustrated in the 
figure by the small plus signs, and the lidar turbulence intensity, 
indicated by the small circles. The interesting result is that the 
turbulence intensities, although scattered, are intermingled, indicating 
general agreement between the 1 idar-measured turbulence intensity and 
the aircraft-measured values. This is particularly true for the lower 
range bins. 

This result is an important observation. It is apparent that 
results from the present study contradict this thinking. It is 
generally thought that the Doppler second moment data will correlate 
with essentially point measured turbulence intensities obtained from 
the aircraft. The fluctuations in the radial wind component time 
history, on the other hand, being values of wind averaged over the 
spatial extent of the range bin, are throught not to necessarily 
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correspond to turbulence measured internal to the volume element. 

As noted earlier, only three range bln values of spectral width 
determined turbulence intensities could be extracted from the lidar 
signal. These values, converted to meters per second in Table 2, range 
from 1.26 to 2.51 m/s; almost a factor of ten larger than values mea- 
sured by the aircraft or computed from tfie lidar data as described 
above. 

In order to investigate the turbulence measurements further, the 
turbulent energy spectra were computed. Turbulence spectra were com- 
puted for each of the eight flight paths and at each corresponding range 
bin, assuming vertical homogeneity. The spectrum computed for each 
range bin for the eight aircraft flights was then segment averaged to 
provide the spectra illustrated by the small plus signs in Figure 15. 
Similarly, spectra for a 2-minute time period begin at the time the 
aircraft enters the range bin, or a region parallel to it, were then 
computed from the lidar data. Note these data are sampled at approxi- 
mately two times per second resulting in a Nyquist frequency of approxi- 
matley 1 Hz, The aircraft data, on the other hand, are sampled at 40 
times per second resulting in a Nyquist frequency of 20 Hz. The spectra 
computed from the lidar data were only five segment averaged. The 
reason for this is illustrated by inspection of Figure 12. At times 
corresponding to some of the later aircraft flights, the radial wind 
measured by the lidar at the higher elevations or higher numbered range 
bins (i.e., approximately Range Bins 16 through 21) were extremely 
intermittent. This is probably due to cloud formation during the later 
runs. Therefore, although these time histories provide a reasonably 
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valid average or mean wind speed- they do not allow a valid spectrum to 
be computed. At the lower range bins (i.e,, Range Bins 11 through 16), 
very good agreement with the aircraft data is observed, Note Range Bins 
9 and 10 were not used because very few aircraft flights descend to that 
height. 

Although the data do not fall on top of one another because of the 
different sampling frequencies involved, the spectra do merge together 
forming a relatively continuous line. This indicated that the distri- 
bution of turbulent intensity in the frequency domain is essentially 
the same for both measurements. The disagreement in spectra at the 
higher range bins is due to increasing noise or decreasing signal -to- 
noise ratio, which is clearly apparent in Figure 12. 

The very good agreement both in turbulence intensity and turbulence 
spectral properties occurring in the clear-air measurements leads to 
the conclusion that computed values of turbulence properties using the 
time history of the lidar-measured winds provide highly meaningful 
results. Although further research is required, this suggests that the 
second moment or spectral width of the Doppler frequency from the lidar 
may not be necessary in order to compute turbulence properties. If this 
is true, the time history of the wind speeds measured by the lidar can 
simply be analyzed for turbulent statistical properties of interest. 


4.0 CONCLUSIONS 


It is concluded that very good agreement between remotely sensed 
winds using a ground-based Doppler lidar and in situ measurements with 
an instrumented aircraft ^s possible. Results show that turbulence 
intensities computed from time histories measured with the aircraft and 
time histories of the radial wind measured with lidar can be analyzed 
statistically to provide turbulence intensities and turbulence spectra 
which agree well with one another. The results further show that the 
second moment data, as presently computed with the NASA/MSFC algorithms, 
do not provide meaningful comparisons with turbulence intensities 
measured with the aircraft. This disagreement, however, must be 
investigated further in terms of the accuracy of the second moment data 
determined by both the lidar hardware and the algorithm for computing 
the second moment. 

Finally, additional insight into making turbulence measurements 
with lidar can be achieved by analyzing the May 10, 1983, results 
when the aircraft data is available from NASA/Langley Research Center. 
Also, an analysis of the signal-to-noise ratio with the object of 
calibrating the turbulence intensity at higher bin numbers with aircraft 
data should be carried out. The result would provide a reliable and 
accurate technique for calibrating the lidar and, thus, of remotely 
sensing turbulence in the atmosphere. 
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ORIGINAL PAGE 5? 
OF POOR QUALITY 



Figure 4 Sample of velocity data along the lidar beam for 150 sequential 
beams (represents 75 sec of data). 
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Figure 7 Temperature variations with height 
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Figure 14 Comparison of lidar-measured winds with aircraft-measured 
winds and computed turbulence (assumed horizontal 
homogeneity); Aircraft * x, Lidar * 0. 
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TABLE 1. VABS Tape Format (ground-based lidar). 
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18-77 Additional groups of three variables for each range gate. 
The number of gates depends upon the pulse width. 


TABLE 2. Values of Turbulence Intensity Versus Integer Index. 


Integer Index 

0 

1 2 

3 

4 

5 

6 7 

8 

Lldar Width (m/s) 

0.5 

0.5 0.63 

0.79 

0.99 

1.26 

1.58 1.99 

2.51 


9 

10 11 

12 

13 

14 

15 



3.15 

3.97 5.0 

6.29 

7.92 

9.98 

12.56 



51 


ORIGINAL PAGE m 
OF POOR QUALITY 


i 


in 

4-» Hi 
C o O O 
•I- «n rt m 

a. 


t 


m k — 

<J o- rv 

Q — o- o 

I — <3 o m 

(/) »oo 


I J I I 
! I ! ' 


e c* 

w> in 


iP c* ci 
m m in 
W !•> 1*1 


; i 
i 

® 43 
**1 O 
a. *n 
o A4 

o o 


tn »4 a} 
0 0*0 

•~i 4> (Tl 

o mm 


o» o* c* o* 

mminifti 

w m n w i 


ONOO 
in h o n : 
mmoc 


i.jjj.jjjjjijjiijjiijiilijjjii 

Air^Airw<-f*A«*»'A*r»A«A»A*A*A*F-r-A-r*f^AJf'-A*r'-A*A»fNr'*A*A*A-r*N 


N N ^ N 
IMJ* O* 
m m m m 
! at m m 


I 


i 

n / mm 
Hnmm 
•v i J> ■£> f*> 
a o i\j 

OOOH 


O O Ol 
mmm 
f*) w w 


i ; 

a- o *n 
*4 *n 4> 
nj ^ o* 
hho 
ooo 


t> C* Ch 
mmm 
n n pi 


o n- c* 
m o* .n 
ooo 
ooo 


a* aj 
o* o» 
4% in 
r» rn 


m *r 

n 43 


O rj 

O C 


o*ooa , c*oooao'0*oj»n*oc?'0 , oioiy 

mmmiruninif'ininifMnininmininin«in«»r 

mannnr|n»nn^nF|f»)nnnwf*]fi«^ 


I 


i 


I l 


I 1 i 


i 
i 

. i 

i i 

N r- ywi j i >r 


l 


K8M»"W)0>»riNffN(l|0(J 
nMjnnmNwn/i-rHOHVHOmo^* - 9 
o — — r\io«H — — oof-m-omc* — — aar-c-o 

NNNonoooi'M r-tno*roo*o — — — i - 

ooNeoono"inNOHNfio o — o o e ■#■ 


o 

at 

v> 

V. 

o 


IW M 


in 

N 


O aj 
5>2 < o 


I 


I t 


I 


£X mo 


nw ONO 
onooN 
c x. in ^ h 

o o ^ r-« n 

O O a C 4 !% 


m o <c B 
AJ «H A* f4 

«oNwm 

<fCHO 


o n- o 
'O^ n 
m *4 aj 
non 
o e o 


l 


•r p 4 
o r* 

43 *r 
43 n 

o o 


I 


I 


I 


I > 


N h 


< 


I 


<0 

a 


<f- 

fl3 

L* 


«♦- 

o 

a> 

i — 

CL 

E 

<o 

tS) 


CO 


CO 

iS 


X3 

a> u 

r— at 

CL VI C 

E <t> 

ID LO QJ 
V) CM «- 
O 
<a • 

•m o 

(O 

a * 

,. u o- <0 

s_ o ’*• <r ' 

«r o 

at at * . . 
BEo«»l 
at & Jj m 
■— i_ — 

co o 

•r- C 10 

•r- 

+-> 

•r- OJ 

■9 E 

I t- r ino 
O l— "t, - - 

Co 0 ^ 


o , 

h^OH 

o , 

* 1 , 0 , 

M l 

, «v 1 in ; 

o , 


« h 

! j 1 1 | ! 

o m ® -4 

-r 4 o <i o o 


H HH ^ 


inNMMMmofr^rmoNwocomno' 
WHFJONO^IMN--»»nV i 4 , KfU^nttO 
cm^ntiOFiy Br\JFiFU)h>innnMHC. 
FJ^OOO«MmO»WiC*rNnjNNOOOhN 

*rrgoHOHOO-c^ir^moHOco^*? 

• • • • ^ 


in a» 
«n <r 

AJ o 

r4 fH 

O 43 


* I » 

' i * 

•r cd m 
m O 45 
43 sO a* 
• ' O O* 
cr n o 


0 #• o m 
o » N <6 
H » in 45 
Q« 4fJ 
O’ + CP 43 


i 

in' a* 

Oj 0 
C ao 
A- C 
4T O' 


I 


I 


o m 4 
Oh : 

I* ■ I 

h ! i 


J 


AJ 


I 


43 

m «& »4 cu q n 44 m 

Bnff^ooNMC'CHNp M4i«f4nninf\jinP w^mnn 
4C^omCi^mhMnHino^mr)4t)OOHmoo«Oir 
inco > ooomrnaQO'OOOOOOc*roo*-io*-ioo43*#**r 


fO 


» ! 1 

*4 m « 
O* 43 m 
N»rs 
a- aj r* 
o 


I 


♦ H4J 

♦ CMP 


I 


O O C" 

o o- <3 
ooo 


O 4? *n d 
c, *4 *4 c* 

•*n m o 

,1-4 IT) O 
AJ H O 


I 

l * < 

O' A* 43 (A *) 

n» «r a/ o c 

HISJCO H 

r-» A- 4* O' o 

4J r4 4> O' O 


«# *• 
A' «A 
<V 


«• O «r m ip 
♦ ho cm 
f in 4 c 1 cj 
* r i*u a . O' 
*- O 43 <C « O* 


rn » 

• M 


o l 

r4 P“ 

i 

o 


t 

H 

1 

1 

■ 

i 

i 

_ * 


I 


o m « 
O N 


n 


I **f 


» fM 

I 


m 

i «** 


I 


I i 


I 


H I 


rl h n- 

i+i#* 


i 


BQOOhirmottiffiOHOoow® 


#.ff OFliW 

5 * *1 

,# i I 

— — »*icc'V><i«>*>*Jcn^a5e3e>-*-* 1 r“rno — fliOOBf-A-ij 1 

— o-r-Cir- — v«\jc* /iHDO-tnrfr-^s>coiM>ONiftiif<. 

oco<f njCOH ooocrc- o-iNOO'!MriPO<i<IO 


I 


< o <y 

I «i 


O' N 

on 1 

1 


if i n o m 

'°n 


I I O 


n 


® o 


fn i! i i i I I H 


n a 

O r-« 


*ronONNOi\jfl^<noi;Fimf')inMn*rin 
nNoOOaOOHDHO «r m r\i ^ mnnoh 
000«r4Hr<l4i 40044 > HH 0 Hp 044 «r 


r4 0 

cn AJ 


• o< H4J *ri*r 

Ol O* © Fl< FI 

,9 ; •; ' 

Ol 


rieecoOHO«f-i'i>n»oo<Hfi'fo^i 
HONH-OHOiic-riO'nmooN^O'OiN 
HDO-fOOOm(-lOO<l»\NOO'NCift<le 


■O 

(O 

e 

o 


in 

X S 

O 


1 *\ 
i fn 
AJ 


I AJ 

I o« 

< m 


I l 


1 


oi u 

W OUJ 

+J sr m 
•r- a ^ 
r UO 
u. <4 

• m ^ 


O 
uj 
» in 


) 


» 


1 ion — 

1 — 

i i tn a- i 

o — 


. . rj 


^4 *H 

1 

1 

| 

fH 

i 

1 

i 

I 

i 

i 

t 


' • J 

: 1 
t 

j||! 


, : : 

! i i . 

! 

! 

r 

i 

i 

1 

. 

in 


i ? 

1 1 . i 

» u i 


!ii‘ 

! | i 
J i 


f- 

UJ 1 

1 


AJ -AJ J 


H N © H 

in *r 


o* ^ n 

in O’ er» m \ m 

O » I AJ 

a 

H 
A- 

, O’ 


A* ! 


I 


I 


•— 21 S— 


1 i • sr z | 

1 -*f x O a 4n in «n| 

I_X:i— »— c>uju-iu-,uj 
m - ui t a * i a 


I 


’ihjo m - — l5lM5l5 0H'Lifln l — 0001^5 
IDVlDUIVlMVIUMIllJUIUIUUO pul UWUJlu- 
la. >Qt li t aODQa, n c- iroocir 


t 


I 


i 


0) ci 
C 'M 

c til 
m 2 r ’ 
si a “ 

• i 


mg (no 

O O Ci) Tf H» OOO 

ooooooooooooooor>oooaoo»4*4 

< < U4 Ui 

QfatKQtofflfaooooifliOQyfrQr«f»:*Qa.Q. 

I ! ! | j I I , ! 1 • ‘ ‘ ’ ' ' 

| 'h ,n»-i-ois«: • - z * I 1 ; ' I 1 1 ■ ’ 1 

iO ! I i uur-a .hi. *- t- 

O j ■ >*~t —• UH «L JkLa 

• i m i" 7 r w >- u h a m * j oh at 

« c . — o. m -v • ■< < « o i— i- 

H- »- JJJJKLLJtlTI l JU 

lHUMHJHJOOUUll.truLHH.HHr 
:: X tlLiui-ia.UUUUJllilll‘UtJJ^-l‘llWUIOUljUliilVLMJU*l «T JJ Ifl v •- H L. m O H Hi 

HHa.a.oio<<<<<o*-i-«<iB«aio,HOOOi>":ooc»cooQPCio<iai.JJii. 

— * or .*1 ^ ^ r.1 — H '"*• ■ " I . ^ 1 — H 

- • Y'H' < :Tr!’r 


or • o c 

*h «r i— h. 

CT) — I LJ Of 

d ar ui « Oft d d 

riiMM«iafJJHpT 
o > y i j o 1 <jiilii:*->-u.i-i 


> m 


aC -t 


’i t 

1 1 i 

1 

O 

o o . ir 

J 


DG 0 ‘ * 


f- T i- * 

V’ 


»IL 


_ I n i«,|,_ i oi . «' r_: * w' - j* _ 

1 • 1 iTi'i : i . "VI- 


I 


I 


i 


/ w» 4 i •» ^ • •* I i 

• > •- a | * 

lit'. t 


i 


Si- 


52 


15t7 


ORIGINAL EH 
OF POOR QUALITY 


£ j g : j J j • j i - j j i j J , ^ ; j J j > j j » j I • j J . j ■ j I • : j j t . j j J > 


0* Ol 

m m in 

m ri 


0 O' 0 

«n m 


«n 4> ** 

O *f M 
O ^ 

*“« 0 AJ 
O VA « 


rsi n 43 
JNC? 
m N O' 


0 0*0 0* 
m tn tn 

w n •h w 

1 1 > 


•00 00000000 * 0000 *( 

im*\*N*rv*\in»nmintfNininmini 

in»*)f»>fn»»»funnnnfinwai 


1 o m 10 

>••••• 

HNhNUI 


<0*0^0 

OMtIH 

n<n s 1 
HO ON 


O O M «3 tf> IA 

■ •••■■ 


ofuo^o «n o 

/ N / h 4 «rON / 41 
f\itr »^330 cr o> N n n 
m m «> c O OO ^nnom 
lA 4 A IA t"» f\J NN it in o h o 


000 

tA tn v> 

mwn 


0*00 
ia m ** 


000 

• • • 

tV tV 


-1*^0 

Of)N 

ICO^ 

intAO 

»TNO 


*r ON 
NO J 3 

O N 0 


4 * Q ^ 43 < 

O *A O r-# I 

O *-> O O 1 

NNC I 
O CHM 


^ « 
<0 O 4> 1 

fn o ^ < 

© CD tA i 


, 4 O- N ON 
nno>m* 4v ^ 1 
, ^ ® N ki H N O I 
I IA O C I 

*AlT»Ai ffVJA^ *T lf( I 


> 4 A O «A 

► ^ 43 

1 as *r m 
I A» *9 *C 
1 O 4* f*- 


NOO 
n fw fu 
OIDP 
43 a* O 
m 43 o 


H40in< 

-4 © m a- < 

OWK* »l 

Hf»- <1 H I 

c *r r- o < 


ON4 A 
NIT N0 
O N O 0 
onmo 
c © r* m* 




• • 4 

«r -* i* 

nniuidO^ m if , h h rt H rt 


tv 



4 ! »f 

H M ^ N ^ O O . 


• » 

1 . 



, i-j n .1 

1 : • ’ 

1 : 1 1 

: ; : ’ ! 



tA fw o v* 

O 43 0 * o 

O MHlf 

vr tv CM o 
o -fln ^ 


— < *r •-* 

C tv 4* 
© a- in 
O **4 ** 
o* r- m 


At f*1 C 
C- 43 *jr 

•r o r 

® N IP 

•* r* O 


J-NlA^COCH^Nr* no 0000000^4: COCIMTC 

£ O S G O © O © ^ U. r* N^OOOOOOOSk’'C OChNN 
*-* <y ^ooOfs'^rv. , OA‘oooocooC'e*vriooor\ o in 
OOtS. 43 C ©OOCOA’ r^mOOCOOOOOOOOOCOC 
f- H 4 v o c 00 ©go 00 ©©©cjooooooogogo 


OMn 
o m cd 
O rX O 

000 

000 


coosc 
O O 43 m O’ 
sT O 43 %r 

O o © ^ -4 
00000 


4 Hfl C 4 HIPN%Nine® 
43*#; H | Hf- A N 4) O 


I I 

i \ i 


: 1 1 " • h . ~~ ’ ■ ! I ! ! 

. i ! 1 ■ 1 1 i ! !V; ! I i ! i ! i : ! 1 

D-rJp!>r)HfDNP^OO^n(fMO^NT 41 a K nN «f 


O‘fHm^flHf 0 N^ 4 OOlftnifti 0 4 N® 4 aNnN«fONPHNHONJ , NHnrJON 
fM^^f,n 4 iOfNiiAr*-OiAO* k riANOtxr^«r 4 :iA 4 >iA»-«o^oo*ro*^Atr»A»Aifw^tA 4 )rv 

ocoa%rff3A»oon^«4-OA-4)*roor-tAi4>r* oco £OonNDNHiniPNOOH*r 00 %r 


n*r Hi*t4 
WNNhr 
^ChNC 


»r r\j 
43 sT 

J ll 


n run 1 
t * 0 1 


O OF*4 1 
43 03 0> 43 i 
O A* «J* O* 


> fp rvi h 1 1 | 

it 1 111 1 , 


1 j- j- 1 

. 11 1 . * 


M H I »* I H 

I f I I 


•■* I •■**•*> 
f II I 


iN •J'NI 
4) a> m < 

N a. O* • 


OnONinC O CO 0 < 4 *T^O**^r 4 CVIf 
1 m n O O N N HNlftiAOO«CPPt< 
»nrtOCx(Mninpioo>r 0 0^1 


H4 OiT 
I 43 O AJ m 

A- 4 : O O 


o ■/ o o n 

4 IP sT 4 ^ 

«r o o o <n 


»■» O N PM 

c? ru sr » c 
0 o o «-• »-«• 


nnrjNH 


1 m rr> t 

I J 


I WWUJ' 
■ ►— •— Or ■ 

i t? n © « 

2 ?UI 

* * o 


> u u U U 'J U U U U C t_> t_> © 

' U' UJ u. UJ W X 111 U. 1 U. U4 U4 Ui UJ u u . 


'UUOU 
I UJ UJ OJ u> 
I l/HP in IP 


J U u u u 

LU li.’ U. 1 u. < 

i/I lA tO tn (O 1 


UUULMJ 
- at ku Uikt 
.v*. in ip '/m/1 


r c c r r r rxrxrrr*:rv 


1 c © u* - 

. a. a* UJ 
► to aj a. . 


1 c pr 1 

' z 1 

3 0 0 1 


in m -♦ 1 

1 z :* » 


■ ! , 1 ■ 

I oc ac ac A- 1 »- 

■ O-* A- 'JJ *-U4 ** ** C I C 

>►- I"*k- I kNinmo .1 O 

ir'u.or-u.or’u.a.a.*-* ** 

1 uiyM^jaiMuiuj nmin to 

’ u <J cr u ^ a o«J*^CQ> a. 

mu t I UT 
I CC O D O O (3 \5 iJ 7 >• X U >■ J 


rrrrrtrrrrrtrr 
• r •- x — t •-» j 


0 

mine 

a. a. •— 

Nin 2 m 

m D«- c. 
a o to *», 
* u 2 W 
> > > -J 


a_ On. 

o o c.* 

X -f «— ~ 

a 1 - “ '/t ^ « 

ju.mii j, 

4 a & o *i 
mu u *’ u 
> > > •> 


7 •-» 
t- X 1 
a a. at 

2 J 4 41 

— T X »- 

— w-i A O. J. 

* 3 iv JJ X f 

♦*V -J <4 4 

-j » j x a 
*»•*>>> 


* lV I £5 • X <» ^ y '» v * 1 


n 


